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Abstract 

We present a QCD-based calculation of the exclusive semileptonic decay — > plDi. 
Using the ideas of Heavy Quark Effective Theory, we discuss the factorization of the 
amplitude. Further, rcsummed Sudakov effects are put in to ensure a consistent per- 
turbative expansion. 
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The study of b-hadron decays has been the subject of considerable interest in recent 
years, as a source of information about CKM matrix elements, fl|] as a laboratory for the 
application of QCD (2|, ||, [|,J6| and the development of theoretical tools such as Heavy Quark 
Effective Theory (HQET)[||. In particular, heavy-to-light decays are interesting because 
they give information on V u b, but they are especially difficult to calculate because of the 
essential presence of strong interactions in the hadronic bound state. Recently, however, in 
[||] a method has been formulated for analyzing the exclusive decay B° — > n l + v\ in the region 
of large hadronic recoil and in the limit of a very heavy b-quark. The approach combines 
HQET, perturbative factorization theorems for exclusive processes ||, and exponentiation 
of Sudakov double-logarithms (7j. In this paper we extend this technology to study the 
exclusive decay Aj, — > plui, also in the limit of large hadronic recoil and large b quark mass. 
In this limit, the method provides an asymptotic regime in which a systematic expansion 
exists with corrections 0{a s {cA.QCD^b)), where c is a calculable constant. For some recent 
model-dependent calculations of decays, see ||. From the phenomenological side, to the 
extent that the very large m& limit is realistic our calculation is motivated by the possibility 
that it could provide another method for extracting V u b from data, complementary to the B- 
meson studies. In addition, comparison of experiment with the predictions of our calculation 
will afford added insight into the applicability of perturbative QCD (PQCD) in processes 
with momentum transfers in the few-GeV range. Further, the analysis of this paper is useful 
for contrasting the decays — > pi v\ and B° — > tt l + ui, which may be helpful in elucidating 
the surprisingly large difference in the A& and B° lifetimes recently seen in experimental 
data |L0|. 

We stress that in this paper we are calculating the purely perturbative contribution to 
the above-mentioned exclusive decay. This would be the dominant one in the limit of a very 
heavy b-quark. However, for realistic values of m& there are could be sizable nonperturbative 
contributions in the form of higher-twist effects . These will be discussed in a future work 
[15]. However, as we shall see below, numerical indications are that a nodified perturbation 
expansion is self-consistent even for realistic values of m&. 

The physical picture of this decay is similar to that for the meson, described in Q. 
Sitting inside the A&, the b-quark decays into a W~ and a fast-moving, nearly on-shell 
u-quark. The u-quark propagates through the remaining hadronic medium, picking up a 
light ud pair over a distance y. Since we are considering an exclusive decay, no gluon 
radiation escapes. If y were large enough (0( A CD )) we would expect considerable gluon 
radiation. Thus, for exclusive processes we expect the outgoing u-quark to propagate only 
small distances (<< 0{ ^ CD ) ) before acquiring a ud pair to form the color singlet which 
will hadronize to form the proton. Then the hard gluons which kick these spectator quarks 
will typically be off-shell by 0{A.QCDfnb)-, large enough that we might expect perturbative 
factorization theorems to apply. 

We will state the arguments and results of our approach and sketch the calculations. 
The calculational details will be published elsewhere (T^]. Our procedure is to judiciously 
combine elements of [|| and 11, [13]. We first identify the sources of long-distance behaviour 
which must be dealt with in order to successfully apply perturbation theory. Soft divergences 
will arise from interactions of soft gluons with the heavy quark, and soft and collinear 
divergences will arise from the interactions of virtual partons travelling in the direction of 
the outgoing proton. As in H, we separate out the soft divergences due to the heavy quark 
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using an eikonal treatment of the b-gluon interactions. All the other divergences we regulate 
using the Sudakov resummation procedure developed in and applied to exclusive decays 
in [11, O]. Of course these do not include the collinear divergences that are collected into 



the universal proton wavefunction. This modification of the perturbative expression takes 
into account the transverse momenta of the partons inside the proton. The resummation of 
the most important contribution behaves like exp[— const x log Q (log ( l ° c ) ) ] in the leading 
logarithmic approximation where b is the conjugate Fourier transform variable to l±. In 
the above Q is a typical momentum transfer which is of the order of m b in our case. 
This Sudakov suppression selects the components of the proton wavefunction with small 
transverse separation, thereby ensuring that the exchanged hard gluon is always offshell 
and consequently perturbation theory for the hard scattering amplitude remains stable 
and reliable. Contributions from additional non-valence Fock states will be suppressed by 
Q( A ^ D ) from standard parton model considerations for hard exclusive processes ||. 

We conclude by calculating the form factor and differential decay rate, combining our 
form factor with model A& and proton wavefunctions. 

We now discuss the factorization of the hadronic form factor. The interaction Hamilto- 
nian relevant to the decay A& — > p I &i is 



H int = 9z Vuh{ u^(i _ 75 )&)(zyxi - 75 )^) (i) 



The amplitude is 



M = %V ub {P{p')\u ltl {l - l5 )b\A b (p)) - 75M) (2) 

V2 „ ' 

We will be concerned with calculating the hadronic matrix element in the limit of 
x ~ 1, where x is defined by 

2 P'P' /Q\ 

x = r (3) 

m b 

where p is the momentum of the A b and p' is that of the outgoing proton, which we take to be 
in the (+) direction. In 0] it was shown that for B° — > ttI + v\ the hadronic matrix element 
could be sensibly written as a convolution of a hard-scattering amplitude with initial- and 
final- state valence hadronic wavefunctions. For the present case the analogous treatment 
is valid. As we discuss below, we can write (suppressing color and spinor indices): 



M,(P,P') = I ^J d " A,(p,p')* Ab (p,k 2 ,k 3 ) (4) 



d A k 2 d 4 k 3 
(2vr) 4 (2vr) 4 

where k 2 and k 3 are light quark initial momenta, ^>\ b is the A b valence wavefunction, and 
Afj, is itself a convolution of the hard-scattering amplitude, H^, and the proton valence 
wavefunction, 'I'p. Written explicitly, 

M»M) = j\d£i] J[d 2 l ± ] J ^^^ P (^p'J ±l )H^ r ,l ±l ,p',p,k 2 ,k 3 )^ Ab (p,k 2 ,k 3 ) 
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[<fc] = <*i-E& Ki« (6) 




[^l^lVil^VlA (7) 



where the £j are the longitudinal momentum fractions of the proton's valence quarks and 
l±i are the corresponding perp components. 

We now briefly outline the steps leading to eqn.[4]. As usual in the treatment of ex- 
clusive processes, to yield a hard-scattering amplitude computable in PQCD the infrared 
divergences must be factored from the hard-scattering amplitude into the wavefunctions. 
The factorization procedure is a straightforward extension of the corresponding analysis in 

We review the salient points of the factorization of soft divergences here. The remaining 
collinear divergences will be controlled by the resummation procedure discussed later. 

Step 1 Decoupling assumption: integrate out the heavy b-quark loop effects, allowing the 
matrix element to be written as 

M, = (P(p')\J,(0)\A b ) CQcD = (P(p')\T(J,(0)V(A°))\A b ) c(q) (8) 

Here V represents a single b-quark line connected to the weak vertex. V is composed 
of a sequence of b-gluon vertices and b-propagators and as such contains the full effect 
of the terms of the Lagrangian Co(b) — igbj ■ Ab, up to loop corrections. It is V which 
now contains the soft divergences. Note that the matrix element is now calculated 
using the light-quark Lagrangian density C(q) only. 

Step 2 Go to the A b rest-frame and introduce the eikonal phase. 

U(A°) = Vexp[-ig / n ■ A{\n^)dX] (9) 

J — oo 

n is a unit vector in the direction of the heavy-quark velocity (here = 5qu)- This 
contains precisely the soft divergences of V(A ) that we wish to remove. The associ- 
ated Feynman rule is shown explicitly in [||. 

Step 3 Define a subracted b-quark line, V-^ m = VU^ 1 , from which the soft divergences 
have been removed. To lowest order this is the difference between the full b-quark 
propagator and the eikonal propagator. 

Step 4 Express the hadronic matrix element in terms of this subtracted propagator, in- 
serting a UU -1 and a sum over a complete set of states to obtain: 

oo 

M, = £ (P(p')\T(J,(0)Vf in (A°) ® n U(A°))\A b ) c(q) (10) 

n=3 

Here the ■ • - ®n ■ ■ ■ denotes a summation over intermediate states with n — 1 light 

partons, integration over their momenta, and summation over other relevant indices. 
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Step 5 Neglect the higher Fock states. What remains after a reorganization of the pertur- 
bative sum is: 

C A^ h A^ h 00 

= J ^4^4 E[W)l r (^(0)^ /in (^°))M<M^(^ ))|A6)] £ (,) (11) 



=3 



where we identify the second factor in the integrand as the A& wavefunction and the 
first factor as from eqn (3). The passage to eqn (4) is justified if the process 
the proton emerges from a short-distance region and the procedure is by now fairly 
standard (see JT^]). 

We next discuss the implementation of the Sudakov suppression a la' Li and Sterman. 



For this purpose, following [11, 12] we work in an axial gauge. The first step is to transform 



eqn (4) into '6' space, which is Fourier-conjugate to Zj 



M ll (p,p')= I [dfc] l[d 2 b] 
where 



r4 1^ ^4 



(2vr) 4 (2vr) 4 
[d 2 b] 



2 , n d 2 b 2 d 2 b 3 



3 J 

(12) 
(13) 



(2vr) 2 (2vr) 2 

[i is the factorization scale, and V p {£i,p', 6j, fj,) are the transformed proton wavefunctions. 
The Vs include all the large logarithmic and double logarithmic radiative corrections at 
large b. These have been resummed and exponentiated |l3| and result in a suppression of 
the distribution amplitude at large b. Explicitly, we write 



Vp(^i,p',bi,n) = exp 



1=1 



dfj, 



y P {t l ,p',w)+0(a 2 {w)) 



(14) 

where w = minj(-J-). The s(^i,bi,mb) are defined in [jl^, |l2| and j q = — — is the axial gauge 

quark anomalous dimension. The bi are the infrared cutoff parameters and we follow the 
"MAX" prescription of |l|], i.e. 



b { = max [61, 62,^3 



(15) 



where b\ = \b 2 — b 3 \. Including the [i dependence of the hard- scattering amplitude we can 
write our final result in the form 



M^(p,p') 



1 



b 2 db 2 



b 3 db 3 



2- 



d,6 



d A k 2 d 4 k 3 
(2vr) 4 (2vr) 4 



(2vr)3 7 1 W J 
*p (& ! P' , w)!!^ (Ci,bi,9,p',p,k 2 ,k 3 ,t ai , t a2 ) * Afc (p,k 2 , k 3 



xexp 



bi, t ai , t. 



a 2 1 



(16) 



where is the angle between b 2 and 63. In the above, t ai ,t a2 are the appropriate scales at 
which the running couplings in the hard scattering amplitude are evaluated and are related 
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to the largest mass scales appearing in H^. The detailed forms are given below. The S are 
generically of the form 



*' dfj, 
i p 



(17) 



where the t\ depend on the hard scattering diagram (see below). 

We now express the wavefunctions in explicit forms more useful for calculations. In the 
above we have neglected any intrinsic l± dependence in the wavefunctions. A more complete 



treatment taking this into account will be given elsewhere [15]. We work in the rest frame 
of the A b with the proton moving off in the (+) direction. Following the notation of other 



authors 16, 17] we write the final-state proton wavefunction: 



^P(^i)a/3 7 ,lmn = ^lmnjfp[(C^f ■ p') a /3(up(p') J5)jV fa, £ 2 , &) 

+{iCa Su p' u} )(up{p'hsl5) 1 T(Ci,a^3)} (18) 
II (^y^ 1 ) (P(p')K(zi)u^(z 2 )d;(z 3 )\0) (19) 



i=i 



where 



2T(£i,fc,&) = 6) + M&>&,&) (2°) 

2^1,6,6) = + <M6,6,&) (21) 

2^(^,6,6) = M&, &,&)- Mb, 6, 6) (22) 

and £1 + £2 + £3 = 1- C 1 is the charge conjugation operator and = | [7^,7^]- /p is 
related to the wavefunction at the origin. It is chosen such that the distribution amplitude 
has normalization ./^[c^i] fipid: = 1- 

We may define the A b wavefunction similarly: 

®A b (fa, h)al3j,lmn = ^Imn [{l ' P + m b )j 5 C] /3ry u(p) a ip(p, k 2 , h) (23) 



d*z 2 e«*-» J d^ 3 e^ z H0\b l a (0)U(A {J ) U f(z 2 )d^(z 3 )\A b (p)) (24) 

Note that our choice of a scalar momentum wavefunction ip(p,k2,k 3 ) (rather than a 
more complicated tensor structure) corresponds to an assumption of decoupling of the spin 
and orbital degrees of freedom in the light-quark system, as discussed in |l8| . 

With our choice of frame we find that to leading order in m b the only nontrivial depen- 
dence of the hard-scattering amplitude on k 2 and k 3 is through their (— ) components. Thus 
it is convenient to define a distribution amplitude $(k 2 , which replaces ip{p, k^ ,k 3 ) and 
in which the (+) and (_L) components have been integrated out: 
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where kf = kf ±kf. It is convenient to change to the dimensionless variables £,77 defined 
W f = Z^rr= and V 



^2 jja 



fc 2 +fc 3 



In terms of these variables we normalize the distribution 



amplitude <f>(^, 77) by 
and define /a 6 analogous to f p : 



d£ / i]drj(f)(^r]) 



af3-y,lmn 



b 



jfA b [(7 • P + mb)l5C} pi u{p)a<j){t,, V) 



(26) 



(27) 



Thus, f\ b is related to the wavefunction at the origin to leading order in a s (m b ) 



The higher order effects can be discussed; see for the meson case. Using the above 
definitions and evaluating the traces, the Dirac structure of the leading terms will be of 
form [up(p')j^(l — j5)u\ b (p)]. The matrix element may then be written (with color indices 
and color wavefunctions suppressed): 



M„ 



1 1 



1 



7) fA b fp[up(p')7^(l -75)«A„(P)] X 



2 (2vr) 3V 4 

rl rl rl poo roo r2ir 

/ / rjdri / d£ / 6 2 (i6 2 / 63*3 / deH^ i ,r ] ,tb 2 ,b 3 ,9,p',p,t ai ,t, 

Jo Jo Jo Jo Jo Jo 



(28) 
(29) 



= [up(p h M (l - !5)uA b (p)]F(x). 

Explicit expressions for He~ for the different diagrams are given below. Note that the 
proton wavefunctions are included in the definition of H. 

We now discuss the calculation of the hard-scattering amplitude, H^. Unlike the B° — > 
1 + V[ case, for — > plvi only a small subset of the full set of 16 tree-level diagrams 
contribute to leading order in m&. In a spacelike axial gauge with the gauge-fixing vector 
having no perp component, the leading diagrams are those shown in figure [1]. After 
contracting all spinor indices with the wavefunction spin-projection operators and all color 
indices to give initial- and final- state color singlets, and after manipulating integration 
variables, their contributions to He~ s are: 



g(al+a2) 
s (a) 

g(bl+b2) 



16(2vr) 4 C> s (^ 2 )as(4) 



-(^(666) + ^(666) + 2T(£ 1 £ 2 6)) x 



64(2ir) 4 CW&M&: 
(6+6) 



K (J^xm 2 b (l-C) V b 3 )-K (J(l-C3)xm 2 b ^b 3 



log 



log?2 



+ log 



log*3 



+ log 



logte 



1 



(^(666) + ^4(666) + 2T(666)) 



^o(V(6 + &)zm27/& 3 )[tf (V&a™^&i) - ^o(V6^r/6i)] 



i=l 



2/3i 



log 



log 1 1 



+ log 



log £3 



+ log 



logt 4 



(30) 
(31) 

(32) 

(33) 
(34) 
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where 



maxfJ (f 2 + &s)xrnfr}, j-] (35) 



t 2 i = max[Jfrxmli]£, — ] (36) 

Oi 



and 



t 3i = max[^/^ 3 xmgr/(l - £), — ] (37) 



1, 



ti = maxfJ (f 2 + CsJsmgTj, »] (38) 

A QCD v 

The remaining ti are defined similarly, as the maximum of b^ 1 and the argument of one of 
the Bessel functions in the associated H. Here j3\ = 33 ^ f , S = log(6Agc£)), and i?o is the 
modified Bessel function. 

The remaining diagrams are either identically zero due to color algebra, as in the case 
of the 3-gluon-vertex diagrams, or are suppressed by 0( C ^ D ) due to Dirac algebra. Note 
in particular that all diagrams involving a subtracted heavy quark propagator (henceforth 
called slashed diagrams) are suppressed by at least 0( C ^ D ) with respect to the leading 
diagrams. It's instructive to contrast this with the meson case to see why. Consider first 
the meson the two tree-level diagrams are shown in figure [2] . The slashed diagram naively 

rvL 



appears suppressed with respect to the other because the heavy quark propagator is : + _ 
0(— ) whereas the light-quark propagator of the other diagram is 7 <yP l +, k } which is naively 
0( A ^ c ). The p 1 in the light-quark numerator is eliminated by the final-state spinors, 
however, leaving ~ O(^) - the same order as the slashed propagator. For the 

baryon, consider the diagrams of figure [3]. Analysis of these diagrams is simplified by the 
fact that their leading terms come from the g^ u pieces of the gluon propagators in this 
choice of axial gauge. Compare the leftmost quark propagators on the top line. Once again 
the slashed propagator appears to be suppressed. This time the leading piece p' of the 
leftmost light-quark propagator in the second diagram is not eliminated, since the adjacent 
light-quark propagator is offshell and shields it from the final-state spinors. Thus our naive 
counting of powers of m& remains correct and the slashed diagram is suppressed with respect 
to the unslashed diagram, in contrast to the meson case. Combining eqns. (22), (23) with 
(24)-(28), we may obtain the expression for the form factor F{x). This is the only dynamical 
difference that we find between the meson and baryon formfactor calculations. 

Next we proceed to the calculation of the differential decay rate. Neglecting lepton and 
proton masses, the differential rate is: 



dT 



dx 



^x 2 (3-2x)|F(x)| 2 |K fc | 2 (39) 



To proceed further we need models for the proton and A& wavef unctions. Several models 
for the proton wavefunction have been proposed in the literature. We adopt the model of 
Chernyak and Zhitnitskii [l(|] (neglecting Q 2 evolution): 



b c P z {i u ^ = lGeV) = 1206e 2 6[11.35(e i 2 +e 2 2 ) + 8.82e 3 2 - 1.68C3-2.94-6.72(e 2 2 -ei 2 )] (40) 
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There is little data available on heavy baryon wavefunctions, but some plausible models 
for A-baryons with unequal mass constituents have been proposed from which we choose 
the following [20]: 



inr 



mr 



2/32(1-7/) 2^(1-0 



(41) 



N is a normalization constant chosen to fulfill eqn.(26) and mi is the light constituent quark 
mass. In the above, (3 is a parameter related to the string tension and is to be fitted to 
data. We present results with two typical values for this parameter. 

From the form of the distribution amplitudes give above and expressions (28), (30), and 
(32), we see that there are no remaining singularities from any endpoint regions. 

Next the differential decay rate may be evaluated numerically for different values of x. 
To show the variation of 4^ with x , we write it in the form 



dr 

dx 



(42) 



R p (x) is plotted against x for x near 1 in figure [4] for two different values of f3, f3 = 0.5 
GeV (solid line) and (3 = 1 GeV (dashed line). In Table I we list the partial decay rates 
for x integrated from x\ to X2- To get an order of magnitude estimate we use V u b ~ 0.003, 
f p = 5.3 x 10~ 3 GeV 2 16], and f\ b ~ f p . The corresponding numerical values for the 
partial rates are given in the last column of Table I. The total width of A;, is approximately 
1.6 x 10~ 12 GeV |2l|| , so that the perturbative partial widths calculated here are of order 
10 -8 of the total A& decay width. This is somewhat smaller than the branching ratios 
for the analogous meson process B° — > 7r/ + f/ calculated in We stress, however, that 
there is considerable uncertainty in the form and parameters of both the A;, and proton 
wavefunctions (see e.g. p2| ) and in the value of f\ b , which could have considerable impact 
on our numerical results. Thus our numerical values for the decay rates should only be 
considered to be rough indicators. 

The branching ratio for A;, — » A c Wi is estimated to be approximately 10 -2 from model 
calculations Q . Accounting for the difference in the CKM elements and scaling down from 
the heavy-heavy decay, we can get an order of magnitude estimate for the branching ratio 
of A;, — > pi 9[ of 10 -8 , which is consistent with our result given the uncertainties. 

We next turn to a discussion of the numerical analysis. We have investigated numerically 
the contributions of different regions in 6-space to determine whether the perturbative 



contribution is dominant. Following the analysis of [12] we look at the contributions to the 
form factor F{x) from different parts of the integration region. This is done by checking to 
see if the dominant contribution comes from the small-6 regions of 6-space. We check this 
numerically by cutting off the each of the 6-integrals at some maximum b c , i.e b%, 62, 63 < b c . 
Our results are shown in figure [5] for j3 = lGeV. We observe that for Kqcd^c > 0.8 the 
value of F(x) is essentially constant, indicating that most of the contribution to the integral 



comes from the region b± , 62 , &3 < 



The curves with the largest values of x appear to 



flatten more quickly than the smaller-x curves at large 6 C , indicating that for more energetic 
outgoing protons the perturbative regime is more dominant than that of less energetic ones, 
as expected. 
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Following [12], we may define b\ji to be the value of b c at which 50% of the total value of 

F(x) has been accumulated. A loose criterion proposed for the consistency of perturbation 

theory that a 2 s (-rr- ) < 1. Our result satisfies this criterion: for x = l,j3 = 0.5 GeV, b^n 

"1/2 1 
is 0.40 and a s = 0.83. Other x and (3 tested give larger values of 61/2, but even the case 

(3 = 1.0 GeV, x = 0.75 gave a s = 0.91. This indicates that most of the value of F(x) is 

accumulated in the short-distance region of the integration space in which we may expect 

the perturbative expansion to be self-consistent. 

In this paper we have calculated the perturbative contribution to the decay A& —* p I v\ 
to leading order in This is the dominant contribution as mj, — ► 00 and we have 

investigated the self-consistency of the perturbative approach for realistic values of m&. 
With our choice of model wavefunctions for the heavy baryon and proton we find that the 
perturbative contribution of A& — » p I vi to the total width of A& is of order 10~ 8 of the 
total width. Our confidence in the self-consistency of the the result is supported by the fact 
that the largest contribution to the form factor comes from the short-distance regions of the 
integration space, as discussed above. We stress again that there is considerable uncertainty 
in the numerical results due to uncertainty about the wavefunctions. However, for at least 
one plausible set of model wavefunctions we obtain values for the decay rates comparable 
to those calculated for B° — > i>i using a perturbative expansion made consistent by the 
introduction resummed Sudakov effects. The issue of the contributions from higher orders 
in ^- (higher twist) is beyond the scope of this paper. 

We are very grateful to George Sterman for many helpful conversations and useful 
suggestions. We would also like to thank Graham Kribs for help with the computer pro- 
gramming. One of us (R.A.) would like to thank the Service Physique Theorique, Saclay, 
for hospitality and support. This work was supported in part by the U.S. Department of 
Energy. 
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Figure Captions 



Figure 1 Leading Tree-Level Diagrams in Spacelike Axial Gauge 

Figure 2 Tree-level Diagrams for B° — ► irl + vi 

Figure 3 A Suppressed Diagram and a Leading Diagram 

Figure 4 R p {x) vs. x for (3 = 1 GeV (dashed line) and /3 = 0.5 GeV (solid line). R(x) 
is in units of GeV -7 

Figure 5 40F(x) vs. b cutoff for = 1 GeV. The curves are x = 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 

intersecting the 6=1 line from top to bottom respectively. b cuto f / is in units of 
l 
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Table 1. 



Process 




x 2 


rS^(GeV) 

(3 = 1 GeV 


rS^(GeV) 
fj = 0.5 GeV 


A b pWi 


0.75 


1.0 


2.4 x 10" 21 


5.7 x 10" 21 




0.85 


1.0 


1.1 x 10" 21 


2.8 x 10" 21 




0.75 


0.90 


1.7 x 10" 21 


4.1 x 10" 21 
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